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' Abstract 

C/^ ■ We discuss the notion of reduction of a special type of explicit solu- 

pH I tions which generalize the solutions appearing in the classical Laplace cas- 

'^I^ ' cade method of integration of hyperbolic equations of the second order in 

, ^ , ■ the plane. We give algorithms of reduction and prove that different natural 

precise definitions of reduction are equivalent. 
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P ■ 1 Introduction 

o ■ 

I Classical methods of integration of linear partial differential equations and systems 

O ■ of such equations (cf. [121 [131 [13) were considerably generalized in the last decades 

^ ■ in [HI [HI [HI [19] . New applications of these methods were found in two-dimensional 

^ . spectral theory [6] and the theory of stochastic systems [H]. The following ex- 

pressions appear as examples of exact solutions in the methods developed in the 
aforementioned publications: 



f{x,y)= aoix,y)ip{x)+ai{x,y)(p'{x) + ... + anix,y)(p^'^\x)+ 
+ bo{x, y)i'{y) + ... + bm{x, y)^^"'Ky) 

with arbitrary functions (p{x), ip{y) and fixed coefficients ai{x,y), bj{x,y). First 
examples of such solutions were found already very early by Euler [8]. For systems 
( [31 [511111 [19] ) and higher-order equations ([l6l[T7lll8j) solutions may appear as sums 
of several expression of the form fll.ll) . Already Darboux and other authors of the 
classical period had noted that the expressions of the form (11.11) sometimes may be 
simplified (i.e. the orders n and m can be made lower) if one introduces new arbitrary 
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functions, for example as 'd{x) = rQ{x)(f{x) + . . . + rk{x)ip^''\x). No general theory 
of simplification (reduction) of the expression of type fll.ip is currently available. In 
this paper we give algorithms of reduction and prove that different natural precise 
definitions of reduction are equivalent. For simplicity we separate the parts of the 
expression (11.11) containing ip{x) and ip{y), and give the following definition: 

Definition 1.1. Euler integral in the plane is an expression of the form 



where v?(x) is an arbitrary function of the variable x and ai{x,y) are given functions 
of the two variables {x,y). 

Below we assume that ai{x,y) belong to some constructive differential field of 
functions in the plane. In order to guarantee correctness of our algorithms and 
results we need to require that we can constructively decide if a given element of this 
field is zero; we also obviously need at least existence of derivatives of the coefficients 
ai{x,y) up to the order needed in the computations. All functions should have a 
common domain of definition — some open subset of M^. The simplest practically 
important example of such a field is the field of rational functions Q{x,y). 

Hereafter we will often call Euler integrals (11.21) simply "integrals" and write 
them as X = L{x, y)(f{x), where 

L{x,y) = ao{x,y) + ai{x,y)D + . . . + an{x,y)D'', D = 

ax 

is a linear ordinary differential operator (LODO). All differential operators in this 
paper will be linear ordinary differential operators in x, i.e. they will include the 
derivatives = d^/dx^ only. Dependence of the coefficients of an operator on the 
both variables x and or on x alone will be explicitly shown after the sign of an 
operator. Everywhere below we will denote linear ordinary differential operators in 
X with upper-case Latin letters and all given functions with lower-case Latin letters; 
lower-case Greek letters will denote arbitrary functions on x. 

Definition 1.2. Euler multiintegral in the plane is an expression of the form 

J = aio{x,y)(pi{x) + an{x,y)ip[{x) + . . . + ain^{x,y)ip^^^\x) + 

+ a2o{x,y)(f2{x) + a2i{x,y)(f2{x) + . . . + a2n2{x,y)ft^\x) + 

+ ako{x,y)ipk{x) + aki{x,y)ip'f^{x) + ... + aknSx,y)iff''\x) = 
= Li{x,y)Lpi{x) + . . . + Lkix,y)(^kix), 

where ipi{x), . . . , ipk{x) are arbitrary functions of the variable x alone and aij{x, y) 
are given functions of the two variables {x,y), 1 <i < k, 1 < j < rij. 

For the sake of brevity we will denote the operator row [Li{x, y), . . . , L^^x, y)) 
as L{x, y), and the column of functions (^(pi{x), . . . , V5fc(a;))* as (p{x): 



I = ao{x, y)'f{x) + ai{x, yW{x) + . . . + a„(x, y)Lp^''\x), 



(1.2) 



J = L{x,y)0{x). 



(1.4) 
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In Petren's thesis [T7j the following linear partial differential operators in the 
plain were considered: 

p— 1 p— 1 

P{x,y) = ^ai{x,y)DlD^ + ^hi{x,y)Dl = A{x,y)D^ + B{x,y)] a^^i ^ 0. 

j=0 1=0 

(1.5) 

Now we will show that any Euler integral 

X = L{x,y)ip{x) (1.6) 

can be an (incomplete) solution of some equation of Petren type P{x,y)X = 0. 

Suppose we have an integral X generated by an operator L{x, y) = ln{x, y)D'^ + 
ln-i{x, y)D'^~^ + . . . + /o(a;, y). We have to find an operator P{x, y) of the form (11.51) . 
such that 

{P{x,y)L{x,y))v{x) = Q (1.7) 

for any ^{x). Since 

P(x, y)L{x, y) = A{x, y)D^L{x, y) + B{x, y)L{x, y) = 
= A{x, y)L{x, y)D^ + A{x, y) [D^, L{x, y)] + B{x, y)L{x, y) = 
= A{x, y) [Ux, y)Dl+^ + /„_i(x, y)D^^ + ... + k{x, y)D,] + 

+A{x, y) [{in)'Ax. y)D^ + {L-iYAx, y)D:-^ + • • • + {kyA^. y)] + 

+B{x, y) [/„(x, y)D^ + Z„_i(x, i/)Dr ' + • • • + h{x, y)] 

and ^{x) is an arbitrary function, then (II. 7p is equivalent to the system 

' A{x,y)ln{x,y) = 0, 
A{x, y) [ln{x, y) + (/„_i);,(a;, y)] + B{x, ?/)/„(x, y) = 0, 

(1-8) 

A{x, y) [Iq{x, y) + (/i)i.(x, y)] + 5(x, y)li{x, y) = 0, 
. A(a;, y){lo)'{x, y) + B{x, y)k{x, y) = 0, 

where li{x, y) are the given coefficients of the operator L{x, y). Since (11.81) is a system 
with (n + 2) linear homogeneous algebraic equations for 2p unknown coefficients 6o, 
6i, . . . , 6p-i, flQ, . . . , ap_i, one easily concludes that (II. 8p always has a nontrivial 
solution if p > and consequently every integral of the form (II. 6p is a solution 
(certainly an incomplete solution in general) of some equation of Petren type (II. 5p 
of sufficiently high degree. 

Analogously for multiintegrals (II. 3p . it is possible to find some equation of Petren 
type satisfied by a given multiintegral. In this case the coefficients of the equa- 
tion (II. 5p satisfy a set of systems of the form ( II. 8p (a system (II. 8p is formed for 
every operator Li{x,y)) with the total number of linear homogeneous algebraic 
equations equal to = J2i=i ^« + where rii is the order of the operator Lj(x, y). 
So if 2p > there is a nontrivial solution of such a system for the coefficients Oj, 6j 
of an equation of Petren type. 
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The linear case considered in this paper can be a basis for a further study of 
exact solutions of nonlinear PDEs integrable by Darboux method: as we know (see 
[21 M HH]) Darboux integrability is equivalent to integrability of the corresponding 
linearized equation by the Laplace cascade method. 

2 Euler integrals 

In order to give the most general precise definition of reducibility of a Euler integral 
(11.21) we indroduce the following natural definition: 

Definition 2.1. The set of all functions Z{I) of the two variables {x,y) which will 
be obtained after substitution of arbitrary (smooth) functions (f{x) into 1^1. S\} will be 
called the function stock generated by some Euler integral I of the form 1^1. 

Any such function stock is an infinite-dimensional linear space; its elements may 
be obviously added and multiplied with constants. 

Definition 2.2. For any Euler integral 1^1. S\} we will call the order n of the leading 
derivative of the arbitrary function f{x) the order of the integral. 

Definition 2.3. We will call a Euler integral / li.ij)) order-reducible if there exists 
another integral 2i = bo{x,y)ip{x) + . . . + bm{x,y)ip^"^\x) of smaller order m, m < n, 
generating the same function stock as the given integral X. 

Note that we do not assume any relation between the functions '^[x) and '(\)[x) 
in X and X\. As we prove below, if an integral is order-reducible then there exists a 
differential relation between '^[x) and il}[x). 

Definition 2.4. We call a Euler integral 1^1.2) operator-reducible if the correspond- 
ing operator L{x,y) may be represented as a composition of operators 

L{x,y) = M{x,y)R{x), (2.1) 

where M{x, y) is another operator of lower order and the coefficients of the operator 
R{x) depend (as shown) only on x. 

In this case the function stocks generated by the integrals I = L{x,y)ip{x) 
and Xi = M{x,y)ijj{x) obviously coincide and there exists the differential relation 
ip{x) = R{x)ip{x). 

Definition 2.5. We will say that a Euler integral ( fi.^j) has no kernel (is kernel- 
irreduciblej, if 2 = L{x,y)(p{x) identically vanishes as a function of the two variables 
only for ip{x) = 0. In the opposite case we will call the integral X kernel-reducible. 

Below we will need the following constructive way of order reduction for any 
Euler integral. 
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Algorithm RI 



The input data of the algorithm is a Euler integral fll.2p . The output of the 
algorithm is a representation of this integral in the form (12.11) with an operator 
M{x, y) of minimal possible order or a proof that the given integral is operator- 
irreducible. 

Let us take two new formal variables yi and ?/2- Consider L{x, yi) and 1/2) as 
two different ordinary operators in with coefficients depending on the parameters 
yi and y2 respectively. One can find the right greatest common divisor of these 
operators rGCD{L{x,yi), L{x,y2)) using the well-known non- commutative version 
of the Euclidean algorithm (cf. for example [Ij), performing divisions of these linear 
ordinary differential operators with remainders: 

L{x,yi) 
L{x,y2) 

Bk-iix,yi,y2) 

So we have rGCD{L{x, yi), L{x, 1/2)) = Bk{x, yi,y2)- 

Case 1: ord-Bfc(x, yi, 2/2) = 0. Then the operators L{x,yi) and L{x,y2) have no 
nontrivial right greatest common divisor and our algorithm terminates. The integral 
I = L{x, y)ip{x) is operator-irreducible. 

Case 2: oid Bk{x,yi,y2) = ordL(x, = ordL(x, 7/2) = n. This means that 
the Euclidean algorithm has terminated on the first division step and L{x,yi) = 
ao{x,yi,y2)L{x,y2), ordao = 0. Consequently L{x,y) = ao{x,y,yo)L{x,yo), where 
yo is some fixed generic point (such that a„(a;, ?/o) ^ and ao{x,y,yo) is correctly 
defined), i.e. L{x, y) is a LODO with coefficients depending on x only, multiplied on 
the left with a function in the two variables (x, y). Thus we have obtained a repre- 
sentation (12.11) with ordM(x,|/) = and R{x) = L{x,yo). As a result we conclude 
that the integral I is operator-reducible to another integral Xi = ao{x,y,yo)ip{x) of 
order zero. 

Case 3: ord i?fc(a;, yi, 2/2) = m, < m < oid L{x,y). Then, if Bf:{x,yi,y2) = 
Bk{x), i.e. does not depend on yi, y2, R{x) = B^^x) so we again obtain the required 
reduction of the integral. In the opposite case we introduce a new formal variable y^ 
and return to the stage (12. 2p of our algorithm RI applying the Euclidean algorithm 
to the operators L(x, 7/3) and Bk{x, yi, y2). Again some operator bI^\x, yi, 2/2, ys) = 
rGCD{L{x,yi), L{x,y2), L{x,y3)) will be found, and in the Cases 1 and 2 (when 
Old Br^^ = or ord-Br'^^ = oidBk) the algorithm terminates; in the Case 3 we return 
to the stage (12. 2p adding a new formal variable 2/4 and finding B^ (x, yi,y2, 2/3, yi) = 
rGCD{L{x, ?/i), . . . , L{x, 1/4)). Since on every cycle for the Case 3 the order of 5^*^ 
decreases at least by one, after a finite number of steps the algorithm terminates. 
As a result the operators M{x, y) and R{x) appearing in the representation (12. ip 
will be found; obviously the order of R{x) will be maximal and the integral Xi = 
M{x,y)ip{x) will be operator-irreducible. 

Note that for an operator-irreducible integral X = L{x,y)(p{x) in the process of 
the work of the algorithm one can find points yi, y2, ■ ■ ■ , ys, s < n + 1, and operators 



= Ao{x, yi,y2)L{x, 2/2) + Bi{x, yi, 2/2), 
= Aiix,yi,y2)Bi{x,yi,y2) + B2{x,yi,y2), 

= Ak{x,yi,y2)Bk{x,yi,y2). 
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Pi{x), . . . , Ps{x) such that 



Pi{x)L{x,yi) + . . . + Ps{x)L{x,y,) = 1. 



(2.3) 



In fact one can take generic i/i, i.e. some points which don't coincide with the poles 
of the coefficients of all operators appearing in the Euclidean algorithm and don't 
make their leading coefficients identically zero. The operators Pi{x) are found using 
the standard reverse substitution in the equations (12.21) . 

One should note that in the general case the number of the points i/i can 
not be made smaller than n + 1. Below we construct an example of such 
second-order integral X = L{x,y)ip{x) that V?/i,?/2, f^GCD{L{x,yi), L{x,y2)) = 
Bk{x,yi,y2), OTdBk{x,yi,y2) = 1, but there exist points 1/2,1/3, such that 
rGCD{L{x,yi),L{x,y2),L{x,y3)) = 5^^^(x, yi, ?/2, Z/s) = 1- 

To this end let us consider the following three-dimensional linear space of second 
degree polynomials: 



V = {p{x) 



a 



1 + P ■ X + 'J ■ x'^l a, /?, 7 G M}. 



Introduce the following two-dimensional subspaces Wy depending on a parameter y 
which are spanned by the basis 



zi=y-l + {y + l)-x, 
Z2 = y ■ X + {y + 1) ■ x'^ 



{Zl,Z2). 



As one can readily check, for any fixed y the functions zi, Z2 are linearly independent 

iWy 

{0}. This allows 



and for every two distinct values of the parameter yi ^ y2, Wy^ and Wy^ have a one- 



dimensional intersection 
us to construct the following LODO 



On the other hand Wy^ fl Wy^ n Wy^ 



L{x,y)if{x) 



z\ z'l 



Z2 Z2 Z2 



z[ 

Z2 Zo 



-1 



2(y + l) 



-D 



2{y + lf 



xy + X + y {xy + x + yY 
L{x,y2)) = Bi{x,yi,y2) = 



with rGCD{L{x,yi),L{x,y2)) = Bi{x,yi,y2) = {yi - y2){{xyi + x + 
yi){xy2 + X + y2)D - {2xyiy2 + 2xyi + 2xy2 + 2x + 2?/i?/2 + l/i + I/2), but 
rGCD{L{x,yi), L{x,y2), L{x,y3)) = 1. So one concludes that in our algorithm 
RI for the constructed operator L{x, y) one can not find two points ?/i, y2 and oper- 
ators Pi{x), P2{x), such that Pi{x)L{x, yi)+P2{x)L{x, 1/2) = 1 although the integral 
X = L{x,y)(f{x) is operator-irreducible. On the other hand for three generic points 
Vi, y2, ys the equality (12. 3p holds. One can easily construct in the same way an 
example on an integral of order n for which one has to choose at least n + 1 points 
yi in order to get (12. Sp . 



Theorem 2.6. The three definitions \2.3l 2.4 and \2. 51 of reducibility of Euler integrals 
are equivalent. 
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Proof. If an integral I is operator-reducible, i.e. if T = L{x, y)ip{x) = 
M{x, y)R{x)(p{x) with some non-trivial R{x) then the integral I has a nontrivial 
kernel which at least contains the space of solutions of the equation R{x)ip{x) = 0. 
Hence X is kernel-reducible. 

If an integral X is kernel- reducible, that is there exists a non-zero function (p{x) 
such that X = L{x,y)(f{x) = then obviously L{x,y) is representable in the form 

L{x,y) = Li{x,y) ( - — ) , 
\ax (fi J 

so the integral X is operator-reducible. 

It is obvious that operator-reducible integrals are order-reducible. In order to 
prove the converse we need the following Lemma: 

Lemma 2.7. Let two Euler integrals X = L{x, y)ip{x) and Xi = M{x, y)il){x) of the 
form M.2\] are given, and Xi is operator-irreducible. Let the following inclusion for 
the function stocks generated by the integrals hold: Z{X) C Z{Xi). Then the order of 
X is greater or equal to the order ofXi, the integral X can be operator-reduced to the 
integral Xi (i.e. L{x,y) = M{x,y)R{x) for some operator R{x)), and the function 
stocks generated by the integrals X and Xi coincide. 

Proof. Applying the algorithm RI to the integral Xi = M{x,y)ilj{x) we find points 
yi, y2, . . . , ys and operators -Pi(x), . . . , Ps{x), such that (12.31) holds for M{x,y), 
since Xi is already operator-irreducible. Substituting the points yi, i = 1, . . . , s, into 
the operator L{x, y) which generates the integral X, we obtain s LODO L{x, yi) with 
coefficients depending on x. 

For any function ip{x) there exists a function ip{x), such that 

L{x,y)(p{x) = M{x,y)'ilj{x), 

since Z{X) C Z{Xi). Hence 

L{x, y,)^{x) = M{x, y.;)'ip{x) (2.4) 

for i = 1, . . . , s. Multiplying the equalities (12.41) with the corresponding Pi{x) and 
adding we obtain 

s s 

^Pi{x)L{x,yi)ip{x) = J^Pi(x)M(x,?/i)V'(x), 

i=l i=l 

or in view of (12. 3p for the operator M{x, y), we obtain 

s 

^Pi{x)L{x,yi)(p{x) = ip{x). 

i=l 

This means that ip{x) is differentially expressible in terms of ip{x): iplx) = R{x)(f{x), 
R{x) = J2iPi{^)^i.^yyi)^ so L{x,y)(p{x) = M{x,y)R{x)(p{x), i.e. X is operator- 
reducible to Xi. Since for every function iIj{x) there exists a solution v^(x) of 
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an ordinary differential equation iIj{x) = R{x)ip{x), we have M{x,y)il){x) = 
M(x,y)R{x)(p{x) = L{x,y)(p{x) for any ipi^x) and the corresponding (p{x). Con- 



The end of the proof of Theorem \2.(A 

Let an integral X = L{x,y)ip{x) be order-reducible to some Xi = Li{x,y)ilj{x). 
We must show that X is operator-reducible. One can suppose Xi to be already 
operator-irreducible since in the opposite case one can use the algorithm RI to find 
an operator-irreducible I2 = M{x,y)C,{x) with the same function stock. Since the 
function stocks generated by Xi and I2 coincide, we can apply Lemma 12.71 to the 
integrals X and I2 and obtain that L{x,y) = M{x,y)R{x), i.e. the integral X is 
operator-reducible. □ 

In view of the obtained result we will call integrals simply "reducible" or "ir- 
reducible" without specifying the precise meanings given in the definitions 12. 3[ 12.41 



Corollary 2.8. The kernel of a Euler integral X = L{x,y)(p{x) coincides with the 
kernel of the operator R{x) , obtained in the algorithm RI. 

Proof. From fl2.ll) one can easily see that Ker R{x) C Kerl. Let (p{x) G Kerl. 
Then = L{x,y)(p{x) = M{x,y)R{x)ip{x) = M{x,y)i/j{x), ip = Rip. From this 
we conclude that ip^x) = 0, since M{x,y) is irreducible. Consequently ip{x) G 
KerR{x). □ 
We are summing up the results about the structure of reducible integrals ob- 
tained so far in the following Proposition: 

Proposition 2.9. Among Euler integrals generating the same function stocks there 
is a unique (up to gauge transformations Lq LQ{x,y)g{x)) irreducible integral 
Xq = LQ{x,y)ilj{x) . All other integrals with the same function stock have the follow- 
ing form: X = LQ{x,y)R{x)(p{x), where R{x) is an arbitrary LODO with coefficients 
depending on x only. 

Proof. Let Z{Xi) = Z{l2). Using the algorithm RI if necessary we may assume that 
Xi and X2 are irreducible. Applying Lemma [2n we see that y) = I/2(x, y)Ri{x), 
L2{x,y) = Li{x,y)R2{x), so Ri{x) = g{x) must be invertible in the ring of LODO 
so it must be a zero-order operator. If X is an arbitrary (reducible) integral and Xq 
is irreducible, we can use Lemma [213 again obtaining the required identity L{x, y) = 



Theorem 2.10. Let X and Ti be Euler integrals of the form M.2\) . Then the func- 
tion stocks generated by these integrals either coincide or have a finite- dimensional 
intersection. 

Proof. We again can assume that the given integrals X and Xi are irreducible. Let 
X = L{x,y)ip{x), Xi = Li{x,y)ip{x) and some function f{x,y) is contained in the 
intersection of the corresponding function stocks. Then the following equality holds: 



sequent ly Z(l) = Z(Xi). 



□ 



and [231 



Lo{x,y)R{x). 



□ 



f{x,y) = L{x,y)ip{x) = Li{x,y)ip{x) 



(2.5) 
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for some ip{x), tp{x). Apply the algorithm RI to the integral I. In the process we 
obtain a finite number of points yi, y2, ■ ■ ■ , Us and operators Pi{x), . . . , Ps{x), such 
that (12.31) holds. Substituting the points yi into (12. 5p . making the appropriate linear 
combination and using (12.31) we obtain 

ip{x) = N{x)i;{x), (2.6) 

where N{x) = '^'^^^ Pi{x)Li{x,yi). Substituting (12. 6p into (12.51) we get 

L(x, y)N{x)ij{x) = L^{x, y)^{x) (2.7) 

or 

(L(x,y)iV(a;)-Li(x,y))^(a;) =0. (2.8) 

If the operator L{x,y)N{x) — Li{x,y) vanishes identically then any function iplx) 
satisfies (12.81) . and, consequently, also (12. 7p . From this we conclude that the function 
stocks generated by X and Xi, coincide (Lemma 12.70 . If the operator L(x, y)N{x) — 
Li{x,y) does not vanish, we deduce from (12. 8p that ip{x) belongs to the kernel of 
the integral, generated by the operator F{x,y) = L{x,y)N{x) — Li{x,y). Applying 
the algorithm RI to F{x,y) one obtains the decomposition F{x,y) = Fi{x,y)S{x), 
where Fi{x,y) is irreducible and ip{x) belongs to the finite-dimensional space of 
solutions of the ordinary differential equation S{x)ip{x) = 0. From (12. 5p we conclude 
that the intersection of the function stocks, generated by the integrals X and Xi is 
the image of this finite-dimensional space under the action of the operator Li{x, y). 
This concludes the proof. □ 
Remark. Obviously the proofs of the Proposition 12.91 and Theorem 12.101 imply 
that we have a simple algorithmic way to check if the function stocks generated 
by two given integrals Xi and X2 coincide. For this we find the corresponding ir- 
reducible operators Li{x,y) and L2{x,y) with the same function stocks using the 
algorithm RI. If the orders of Li{x,y), L2{x,y) do not coincide then the function 
stocks Xi and X2 are different. If the orders of Li{x,y) and L2{x,y) coincide, we 
divide Li{x,y) by L2{x,y) on the left: 

Li{x,y) = L2{x,y)g{x,y) + B{x,y), oidB < ordLi = ordL2. 

The function stocks of Xi and X2 coincide if and only if B{x, y) = 0, gy = 0. In order 
to find the dimension of the intersection Z(Ii) n ^(X2) we repeat the procedure 
described in the proof of Theorem I2.im i.e. we form the operators in the left hand 
side of (12. 8p and apply the algorithm RI to split off the maximal right divisor S{x) 
of the operator F{x,y). The dimension of the intersection Z{Ii) fl Z{l2) will be 
equal to the order of the operator S{x). 

3 Multiintegrals 

In this Section we prove the equivalence of two possible definitions of reducibility of 
Euler multiintegrals (II. 3p and give the algorithms of reduction and other necessary 
technical stuff. 
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Definition 3.1. The function stock Z{J'), generated by a multiintegral M.^\] is the 
set of all functions of two variables (x, y) which is obtained after the substitution of 
arbitrary (smooth) functions ^pi{x), . . . , (pk{x) into lil.3\) . 

Definition 3.2. We will call a multiintegral U.3\) operator-reducible, if one can 
represent its operator row [Li{x, y), . . . , Lk{x, y)) as a composition of operator ma- 
trices 

I Ru{x) ... Rik{x) \ 
{Li{x,y),...,Lk{x,y)) = {Mi{x,y),...,Mp{x,y)) : ■.. : 

y -Rpi(x) . . . Rpk{x) J 

(3.1) 

where Rij{x) are differential operators with coefficients depending on x only, Mi{x, y) 
are differential operators with coefficients depending on {x,y); the matrix [Rij{x)) 
is required to be non-invertible and p < k. 

Below we will briefly represent the formula (13.11) as 

L{x,y) = M{x,y)R{x). (3.2) 

Here and below R{x) stands for a matrix of differential operators (not necessarily a 
square matrix). Invertibility of such a matrix is understood in the algebraic sense, 
i.e. as existence of another operator matrix S{x), such that RS = SR = 

Definition 3.3. We will call a multiintegral 1^1. 3\) kernel-irreducible (or we say that 
the multiintegral has no kernel) if J = Li{x,y)ipi{x) + . . . + Lk{x,y)ipk{x) vanishes 
as a function of two variables only for ipi{x) = 0, i = 1, . . . , k. In the contrary case 
the multiintegral J will be called kernel-reducible. 

Note that we do not find appropriate to introduce here an analogue of order- 
reducibility of multiintegrals. As we will see below, even irreducible multiintegrals 
with the same function stocks admit a representation (11.31) with operators Li{x,y) 
of arbitrary high order. This follows from the fact that for arbitrary multiintegral 
one can find representations (13. 2p with invertible matrices R{x) of operators Rij{x) 
of arbitrary high order. 

One can perform the following operations on the matrix R{x): transposition of 
two rows and addition to one of the row of another row multiplied on the left with 
an ordinary differential operator with coefficients depending only on x. 

The following algorithm uses these operations for reduction of the matrix R{x) 
to a special convenient form. 

Algorithm RM 

First of all, note that the matrix R(x) in the representation (13. 2p may be con- 
sidered to be free from zero rows and zero columns. 

If the first column has only one non-zero element, transpose the rows and put 
it into the first (uppermost) place, then start processing of the second column. If 
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there are several non-zero elements in the first column, we perform the division of 
one of them with another one with the remainder, for example 



Rn{x) =T{x)R2i{x) + Q{x) 



Then, if one subtracts the second row of the matrix R{x) multiplied on the left 
with the differential operator T{x) from the first row, one obtains an equivalent 
matrix, which has the operator Q{x) of lower order than the previous entry Rii{x). 
It is obvious that such operations correspond to the operations of the Euclidean 
algorithm (12. 2p and result in an equivalent matrix with R^^^{x) = rGCD{Rii, R21), 
and r'^i = 0. Consecutively applying this procedure to the other non-zero elements 
of the first column we obtain a matrix which has only one non-zero element R^y^ (x) = 
rGCD{Rii, . . . , Rpi) in the first column. 

Next we proceed to the second column. If it has some of the entries R22, R32, • • • , 
Rp2 different from zero, we apply to them the same procedure as above and obtain 
= rGCD{R22, . . . , Rp2). If, contrarily, all of R22, • • • , i?p2 vanish, we start 
processing of the elements R23, R33, ■ ■ ■ , Rp3 of the third column. Obviously this 
procedure is a direct non-commutative analogue of the standard Gauss elimination 
algorithm, so in a finite number of steps we obtain a matrix of the form 



Rst{x) 



( Ru{x) 















* * 

R3q3ix) 





V 







which will be called below echelon matrix. In a more general case (which is of no 
interest to us) echelon matrices may have a few zero starting columns. First non-zero 
elements of each row Rii{x), R2q2{x), R^q^i^x), . . . will be called (as usual) pivots. 
The number of non-zero rows of the obtained echelon matrix is called the rank of 
the initial operator matrix. As proved in P, 0] , the rank is an invariant of a matrix 
and does not depend on the method of its reduction to an echelon form. 

The procedure of reduction to the echelon form is equivalent to multiplication 
of the initial matrix R{x) on the left with some invertible matrix T(x). Therefore 
in the formula (13. 2p we can substitute M{x,y) with M{x,y)T~^{x) and R{x) with 
Rst{x) = T{x)R{x)] so we can assume below that the matrix R{x) already is echelon 
matrix. 

It is easy to see that the matrix R{x) is invertible if and only if it is a square 
matrix and after its reduction to echelon form we obtain an upper-triangular matrix 
with non-zero diagonal elements which are invertible elements in the ring of differ- 
ential operators, that is if they are zero-order operators (functions). The elements 
strictly above the main diagonal may be arbitrary differential operators. 

Analogously to the algorithm RI we give below another algorithm which extracts 
the "maximal" right factor R{x) from a given multiintegral, i.e. for a given operator 
row L{x, y) a representation (13.21) will be found, with a row M{x, y) defining an 
operator-irreducible multiintegral with the same function stock. 
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Algorithm RMI 



The initial data of the algorithm is an operator row {Li{x, y), . . . , Lk{x, y)). 
According to the decomposition (13. ip with an echelon matrix R{x) we see that 
one has to find operators Mi{x,y), Rij{x), such that 



Li{x,y) = Mi{x,y)Ru{x), 

L2{x,y) = Mi{x,y)Ri2{x) + M2{x,y)R22{x), 



(3.3) 
(3.4) 



Using the algorithm RI we find the decomposition (13.31) . where Rii{x) is a LODO 
of maximal order and Mi{x,y) is irreducible. Also as a byproduct of the algorithm 
RI we get points yi, . . . , yg and operators Pi{x), . . . , Ps{x), such that 



Pi(x)Mi(x, yi) + . . . + P,{x)M,{x, ys 



(3.5) 



Now we should find Ruix), M2{x,y), R22{x) such that the equality (13.41) holds. We 
obtain this decomposition in three stages. 
Stage A. Finding R22{x). 
Substituting the found points yi into (13.41) we obtain 

L2{x,yi) = Mi{x,yi)Ri2{x) + M2{x,yi)R22{x). 

Multiplying these identities with the operators Pi{x) on the left and adding, we get 



^Pi(x)Mi(a;,2/i) 
.«=i 

or, in view of (13.51) . 



Ri2{x) + 



^P^(x)M2l 



x,yi 



i=l 



R22{x) = y^^Pi{x)L2{x,yi 



i=l 



(3.6) 



Rl2{x) + 



^Pi{x)M2{x,yi 



i=l 



R22{X) = J2 

Pi{x)L2{x,yi). (3.7) 



1=1 



Multiplying (13.71) on the left with Mi{x,y) and subtracting from (13.41) . we arrive at: 

[M2{x, y) - Mi(x, y) J2i P^{x)M2{x, y^)] R22{x) 

= L2{x, y) - Mi{x, y) Y,i Pi{x)L2{x, yi). 



(3.8) 



The right hand side of this equality is some known operator N{x,y). Hence in 
order to find R22{x) one has to apply the algorithm RI to N{x,y), obtaining the 
decomposition 

K{x,y)R22{x) = N{x,y) (3.9) 

with R22{x) of maximal order; the obtained operator K{x,y) is irreducible. 
Stage B. Finding Ru^x). 
Rewrite (13. 7p as 



i?i2(x) + S{x)R22ix) = T{x) = ^Pi(x)L2(x,?/i) 
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(3.10) 



where the operators T{x) and R22{x) are known. The operator S{x) can be arbitrar- 
ily chosen (this corresponds to multiphcation of the echelon matrix R{x) on the left 
with an invertible upper-triangular matrix) and define the corresponding Ri2{x). 
Stage C. Finding M2{x,y). 
Recall that (13.81) was obtained from (13.71) and (13. 4p . Now we do the inverse 
operation: multiply (l3.1Up with Mi{x,y) on the left and add (13.91) to it: 



Ml {x, y)Ru (x) + [Ml (x, y)Six) + K{x,y)] R22 (x) = L2 (x, y ) . 

Thus we obtain the required equality (13. 4p with M2(x, y) = Mi(x, y)S{x) + K{x, y). 

Note that since the obtained K{x,y) is irreducible we can find points j/2,1, 1/2,2, 
. . . , ?/2,s2 and operators P2^i{x), . . . , P2,s2{x), such that 

Y,P2A^)K{x,y2,^) = l. (3.11) 

i=l 

If the number of operators Li{x, y) in the row L{x, y) is greater than two, we can 
analogously find the required operators Mj{x,y), Rij{x), j > 3. We demonstrate 
this below for j = 3. 

In the equality 

L^ix, y) = Mi{x, y)Ri^{x) + M2(x, y)R2^{x) + M3(x, y)Rz^{x) (3.12) 

the operators Mi(a;, y) and M2(a;, y) are known. Using the points ?/i, . . . , ys chosen 
above and (13.51) we arrive at the following analogues of (13. 7p and (13. 8p : 



K{x,y)R2z{x) 



M3(x, y) - Mi(x, y) P^{x)M^{x, yi) 

i 

= Lsix, y) - Mi(x, y) ^ Pi{x)LAx, yi) 



R33{x) = (3.13) 



Using now the points y2i, ■ ■ ■ , y2s2 and the equality (13.110 . we can cancel the term 
K{x,y)R23{x) and obtain an analogue of (13. 9p : 

K3{x,y)Rs3{x) = N3{x,y) 

with the known operator N^lx, y). This allows us to find i?33(x) using the algorithm 
RI. Repeating stage B we find Ri^^x), R2z{x) with some freedom, analogous to the 
freedom of choice of S{x) in (I3.10p . Repeating the stage C one finds M^^x.y). 

Remark. It may happen that the right hand side of (13.80 vanishes. Then we 
set R22{x) = and find -Ri2(a;) from (13.70 . Thus M2{x, y) is actually absent in (13. 4p 
and we can set -R33(x) = on the following step in (13.120 : in this case the pivot 
element becomes R23{x). Thus Rij{x) becomes a non-square echelon matrix. The 
number of operators Mi{x,y) in this case will be smaller that the number of the 
initial operators Li{x,y), i.e. p < k. 

As the result of the algorithm RMI we obtain an operator row M{x,y) = 
{^Mi{x , y) , . . . , Mp{x , y)) and an operator echelon matrix R{x) in the representa- 
tion (13. ip . It is not obvious that the complete multiintegral J7i = M{x^y)^{x) 
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is operator- irreducible. This fact will be proved below (Corollary 13. 9p . Inclusion 
Z{J) C Z{Ji) is obvious. The inverse inclusion follows from the simple fact that 

the echelon matrix R{x) obtained in the algorithm RMI has no zero rows so the 
system of linear ordinary differential equations R{x)(p = ijj for the unknown ip is 
solvable for any right hand side tp. 

Remark. If the given multiintegral was operator-irreducible then the algorithm 
RMI will given the decomposition (13. ip with k = p and a triangular invertible 
matrix R{x) as the output. The ambiguity of the choice of the operator S{x) 
on stage B will correspond to the possibility to write a representation (13. ip with 
arbitrary triangular invertible matrix R{x) for any multiintegral JT. Note that the 
requirement of non-invertibility of the matrix R{x) without the requirement p < k 
in the definition 13.21 is not sufficient: it is easy to give an example of decomposition 
(13. ip with p > k (so with a non-invertible matrix R{x)) for any multiintegral = 
L{x,y)(p{x). 

Proposition 3.4. The multiintegral (Mi(x, y), . . . , Mp{x, y)) obtained in the process 
of the work of the algorithm RMI is kernel-irreducible. 

Proof. Let (^ipii^x), . . . , ipp^x)^ be a nonzero element of the kernel of the multiintegral 
M = {Mi{x,y),...,Mp{x,y)): _ 

Mip = 0. (3.14) 

By construction of the operators Mi{x, y) they satisfy 



L = MR. 



(3.15) 



First we carry out the proof of the Proposition for the number of elements in the 
operator rows L and M equal to 2. In this case (I3.14p means that Mi{x,y), ipi{x) , 
i = 1,2, satisfy the equation 



= M^{x,y)^l + M2{x,y)iJ2. 
The operators Mi{x,y), M2{x,y) satisfy in turn 

L2{x,y) = Mi{x,y)Ri2{x) + M2{x,y)R22{x) 



(3.16) 



(3.17) 



where Ri2{x), R22{x) are found in the algorithm RMI. Now we carry out with the 
equality (13.160 the same operations as we did in the algorithm RMI for the equation 
(I3.17P in order to find R^i^x), R22{x), i.e. we take the same points yi and the same 
operators Pi{x) which were found in the algorithm RMI, then we compose the same 
linear combination of the equations that in (13. 7p . We get 



ijji{x) + 

From (I3.18P and (13.160 we have 



5^P,(x)M2(x,; 



M^) = 0. 



(3.18) 



M2{x, y) - M,{x, y) ^ P^{x)M2{x, y,) 



(3.19) 
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(an analogue of fl3.8l) in the algorithm RMI). Thus 'ip2{x) = 0, since the algorithm 
RI gave an irreducible integral K{x, y) in the left hand side of this equality. 

Then from f l3.16p we get Mi{x,y)ipi = 0. But Mi{x,y) is irreducible by con- 
struction so ipii^x) = 0. 

Thereby the Proposition is proved for the multiintegral = 
{M,ix,y),M2ix,y)). 

In order to prove the Proposition for the general case J7 = 
(Mi(a;, y), . . . , Mp{x, y)) we analogously consider the equality 

= M^{x,y)Mx) + ■■■ + Mp{x,y)i)p{x), (3.20) 

where ('?/'i(a;), . . . , il)p{x)) is an element of the kernel of the multiintegral J instead 
of the equality 

Lk{x,y) = Mi{x,y)Rik{x) + . . . + Mp{x,y)Rpk{x). (3.21) 

Applying to (13.201) all operations described in the algorithm RMI one arrives at 
Kp{x,y)ipp{x) = 0, where Kp{x,y) is irreducible, which implies i'pi^x) = 0. Per- 
forming the reverse run of the algorithm we find consecutively ipp-i{x) = 0, . . . , 
V^i(x) = 0. □ 

Corollary 3.5. The kernel of the multiintegral J = (Li(x, y), . . . ,Lfc(x, y)) coin- 
cides with the kernel of the matrix R(x) obtained in the algorithm RMI. 

Proof. The algorithm gives us L{x,y)(p{x) = M{x,y)R{x)(p{x), where M{x,y) de- 
fines a kernel-irreducible multiintegral. If <f{x) belongs to the kernel of the ma- 
trix R{x) then (p{x) obviously belongs to the kernel of the multiintegral J7. Con- 
versely if (p{x) belongs to the kernel of the multiintegral, i.e. L{x, y)0{x) = 0, then 
M{x, y)i/j{x) = 0, ip = R{x)(p{x) and in virtue of the kernel-irreducibility of M{x, y), 
ifj = R{x)0{x) = 0, so ip{x) belongs to the kernel of the matrix R{x). □ 

Proposition 3.6. A multiintegral J has an infinite- dimensional kernel if and only 
if the matrix R{x) obtained in the algorithm RMI is non-square (the number of its 
columns is greater than the number of its rows). 

Proof. Suppose that the matrix R{x) in the decomposition (13. 2p obtained by the 
algorithm RMI is a square matrix. By Corollary 13.51 the kernel of the multiintegral 
coincides with the kernel of R. The dimension of the space of columns {ipi{x), 
ipk{x)) which are solutions of the matrix differential equation R{x)(f{x) = 
is finite. In fact, since R{x) is a square echelon matrix without zero rows, we see 
that ipk{x) satisfies the equation Rkk{x)ipk{x) = and consequently belongs to a 
finite-dimensional space. Substituting any of the found (pk{x) into the previous 
equation Rk-i,k-i{x)(pk-i{x) + Rk-i,k{x)ipk{x) = 0, we find that (fk-i also belongs 
to a finite-dimensional space. Then analogously we define all the other fk-2{x), 
. . . ,ipi{x), each of them belong to a finite-dimensional space of solutions of some 
linear ordinary differential equation. Therefore the dimension of the kernel of the 
multiintegral is finite and equals the sum of the orders of the differential operators 
on the principal diagonal of the matrix R. 
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Suppose now that the matrix R(x) is non-square. In any case the algorithm RMI 
outputs an echelon matrix without zero rows and, if it is non-square, the number 
of its rows is smaller than the number of its columns. We show that the kernel of 
such a matrix is infinite-dimensional. Take the first row of R{x) with a non-diagonal 
pivot and consider the previous row, which has the first non-zero elements Rqq ^ 0, 

g+i, . . . , Rq,q+r SO that Rg+i^q+r+i 7^ is the pivot of the next row, r > 1. Choose 
the following functions fi{x), . . . ,ipk{x): set ^i{x) = 0, if i > g-fr; the other ^Pq[x), 
. . . ,(pq-^.r{x) should bc taken from the space of solutions of the differential equation 
Rqq(Pq{x) + . . . + Rq^q+r'^q+r{x) = 0. Thls spacc Is obviously infinite-dimensioual since 
the equation contains at least two unknown functions (pi{x). Substituting the chosen 
iPq{x), . . . ,ipq^r{x) luto the prcvious equations of the system R{x)ip{x) = 0, we find 
the other ipq_i{x), ...,ipi{x). So we have found an infinite-dimensional subspace 
of the solution space of the system R{x)ip{x) = 0. The functional dimension (the 
number of free functions of one variables other than constants) of the complete 
kernel, as one can easily conclude from the considerations above, is equal to the 
difference of the numbers of the columns and rows of the matrix R{x). □ 

Lemma 3.7. Let an operator matrix R{x) of the size k x p with k < p be given. 
Then R{x) is non-invertible if and only if it has a nontrivial kernel. 

Proof. It is easy to see that nontriviality of the kernel implies that R{x) is non- 
invertible. 

Conversely, let R{x) be non-invertible. Reduce R{x) to echelon form. Multi- 
plying the rows of the resulting R{x) with appropriate non-zero functions we can 
without limitation of generality assume that its pivots of order zero are equal to 1. 

Case A: let k = p, i.e. the matrix R{x) is square. Then among its pivots Ru{x), 
. . . , Rpp{x) at least one is not equal to 1 (otherwise R{x) is invertible). Take the first 
pivot different from 1. Let it be Rmm{x). The differential equation Rmm{x)(p{x) = 
always has a non-zero solution <f{x). Consequently the kernel R{x) contains a non- 
zero element ip{x) = {{pi{x), . . . ,ipp{x)y, where 

(Pm+l{x) = ... = ipp{x) = 0, 
iPmix) = ^(x), 

iPm-l{x) = -Rm-l,m{x)'^m{x), 

<^m-2{x) = ~Rm~2,m-l{x)(Pm^l{x) - Rm-2,m{x)'^m{x) , 

case B: now we assume R{x) to be non-square, i.e. k < p. Take the first row of 
R{x) with the pivot Rmm such that the next pivot is not diagonal (or take the last 
row). We have the entries Rmm 7^ 0, Rm,m+i, ■ ■ ■ , Rm,m+r of this row, such that 
Rm+i,m+r+i 7^ is the next pivot, r > 1 (or m = p so we take all elements of the 
last row) . If the pivot Rmm is not equal to 1 , we can easily find a nontrivial element 
of the kernel R{x), using the argumentation of the case A. If we have r + 1 > 2 then 
we can find an infinite-dimensional kernel of R{x) using the argumentation of the 
proof of Proposition 13.61 □ 
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Theorem 3.8. A multiintegral J is operator-reducible if and only if J is kernel- 
reducible. 

Proof. Applying to = L{x,y)if{x) the algorithm RMI, we obtain the represen- 
tation L{x,y) = M{x,y)R{x). If J' is kernel-reducible, then according to Corol- 
lary 13.51 the matrix R{x) has a nontrivial kernel and is non-invertible. So is 
operator-reducible . 

Conversely, let be operator-reducible, i.e. we have the corresponding decom- 
position L{x,y) = M{x,y)R{x), where R{x) is a non-invertible matrix, such that 
the number of its rows is smaller or equal than the number of its columns. Ac- 
cording to Lemma 13.71 the operator matrix R{x) has a nontrivial kernel. So is 
kernel-reducible. □ 

From this Theorem and Proposition 13.41 we immediately get 

Corollary 3.9. The obtained in the algorithm RMI multiintegral J = M{x, y)'ip{x) 
is operator-reducible. 

Theorem 3.10. Let two multiintegrals J\ = [Mi{x,y), . . . , Mp(x,y)) , J2 = 
(^Li{x,y), . . . , Lk{x,y)) be given, and Ji is operator-irreducible. Let the function 
stock generated by the multiintegral J2 be contained (as a subset) in the function 
stock generated by J\. Then 

{Li{x,y),...,Lk{x,y)) = {Mi{x,y), . . . , Mp{x,y))R{x) (3.22) 

with some matrix R{x) of operators whose coefficients depend only on x. Moreover 
the rank of the matrix R{x) coincides with the number p of its rows if and only if 
the function stocks generated by J\ and J2, coincide. 

Proof Since we have Z{J2) C Z{Ji), then for every column (^(x) = {ipi{x), 
. . . ,ipk{x)y one can find another column ip^x) = {ipi{x), . . . ,ipp{x)Y, such that 

(Mi(x, y), . . . , Mp(x, y))ij = y), . . . , y))^. (3.23) 

Reproduce the procedures of the algorithm RMI for the equation f l3.23p with the 
obvious modifications. Since Mi(x, ?/) is irreducible, one can choose the values 
. . . , ys and operators Pi{x), . . . , Ps{x) so that 

Piix)Mi{x,yi) + . . . + Psix)M,{x,ys) = 1. 
Therefore from fl3.23p we can obtain 

s s 

V^i(x) + J2Pii^)M2{x,yi)Mx) + ■■■ = 5^Pi(x)Li(x,y,)(^i(x) + . . . (3.24) 

i=l 1=1 

Applying the operator Mi (x, y) to the both sides of this equation and subtracting 
the result from fl3.23p we arrive at an analogue of the equality fl3.13p of the algorithm 
RMI. Proceeding with the algorithm we finally obtain an equality of the form 

k 

Q{x,y)^p{x) = ^Tj{x)ipj{x). (3.25) 
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Since by the assumption of the Theorem Jx is operator-irreducible, the obtained 
operator Q{x, y) is also operator-irreducible so choosing some appropriate points yi, 
■ ■ ■ , ijt and operators Pi(x), . . . , Pt{x) one has 

Pi{x)Qi{x,yi) + . . . + Pt{x)Qi{x,yt) = 1. 

Then from (13.250 we can easily deduce 

k 

^p(^) = ^Rpj{x)'Pj{x)- (3-26) 
i=i 

The found operators Rpj{x) give us precisely the last row of the matrix R{x). Its 
previous rows can be readily found from the obtained in the process of the algorithm 
equalities of the form fl3.24p . Substituting fl3.26l) and the analogous expressions 
for ipp_i{x), . . . , ipi{x) into (13.231) we obtain the required formula (13.221) since the 
functions (pi{x) were chosen to be arbitrary. 

From (I3.22P one concludes that the function stocks generated by the multiinte- 
grals in question coincide iff the system of linear ordinary differential equations 

R{x)ip{x) = ip{x) (3.27) 

is solvable for every column 'ip{x). Reduce R to echelon form (this is equivalent to 
multiplication of the both sides of fl3.27p on the left with some invertible operator 
matrix S{x)). The system S{x)R{x)ip{x) = S{x)4'{x) is solvable for any ip{x) if 
and only if the echelon matrix S{x)R{x) has no zero rows. This proves the last 
statement of the Theorem. □ 

Note that contrary to the case of Euler integrals, a strict inclusion Z{J'2) C Z{J'i) 
is possible for multiintegrals. The simplest case can be given by any multiintegral 
of the form J7i = [Li{x,y), L2{x,y)){ipi{x), (p2ix)y where both Li(x,y), L2{x,y) 
generate irreducible integrals with different nonintersecting function stocks Z{Liipi), 
Z{L2(f2)- First of all we deduce that J7i is kernel-irreducible and its function stock 
coincides with the direct sum Z{Liipi) © Z{L2^2)- Take now J2 = Li{x,y)ipi{x). 
We see that Z{J2) = Z{L2f2) is strictly contained in Z{J'i). 

This example motivates the introduction of the following notion 

Definition 3.11. // an operator-irreducible multiintegral J = 
(Li(x, y), . . . , Lfc(x, y)) is given, then its submultiintegral is any multiintegral 
J of the form 

J = {Li{x,y),. . Lm{x, y))^{x) = {Li{x, y),...,Lk{x, y))N{x)ip{x), 

where N{x) is an arbitrary matrix of LODO with coefficients depending on x. As 
obvious from the results proved above, this is equivalent to the requirement Z{J') C 

ziJ). 

Remark. As we can see from the proof of Theorem l3.10l there exists an algorith- 
mic way to decide for any two given multiintegrals J7i and 1/2, wether the function 
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stocks generated by them coincide or are different, or wether we have a strict in- 
clusion of one stock into another and what is the dimension of their intersection 

To this end we find using the algorithm RMI the irreducible represen- 
tations of the given multiintegrals Ji = (^Li{x,y), . . . , Lk{x,y))(f{x), J 2 = 
[Mi{x, y), . . . , Mp{x, y))'tp{x) (with the same function stocks). Write the following 
formal equality 

{M^{x,y),...,Mp{x,y))ij{x) = {Li{x,y), . . . , Lk{x,y))^{x), (3.28) 

where ip{x) = {iIji{x), . . . ,iljp{x)y, ^p{x) = {(pi{x), . . . ,(pk{x)y. Applying the proce- 
dures of the algorithm RMI to the equality (13.281) . as in the proof of Theorem 13. 101 
we find a matrix R{x) such that 

^{x) = R{x)ip{x). (3.29) 

From (I3.28P and (13.291) one deduce that the intersection of the function stocks of J7i 
and J'2 coincides with the set L{x, y)0{x), where (p{x) are solutions of the equations 

M{x,y)R{x)(p{x) = L{x,y)(p{x). 

Thus if we have the operator identity 

L{x,y) = M{x,y)R{x), (3.30) 

then either Z^J^) C Z{J^ so by definition J7i is a proper submultiintegral of J2 (in 
the case when the rank of R is smaller than the number of its rows p), or Z{J'2) = 
Z{J'i) (if the rank of R equals the number of its rows p). If the equality (13.301) does 
not hold we form a new multiintegral jTs = K{x, y)^{x) with K{x, y) = L{x, y) — 
M{x,y)R{x). Applying to jTs the algorithm RMI, we obtain the representation 

K{x,y) = N{x,y)Ri{x). (3.31) 

The intersection of the function stocks of J7i and J'2 coincides with the set 
{L{x,y)^{x) I ^(x) e Keri?i(a;)} so by Proposition 13.61 it is finite-dimensional if and 
only if the matrix Ri{x) is a square matrix. In the contrary case (when the num- 
ber of its columns k is greater than the number of its rows pi) Zi^Jq) and Zi^J^ 
has an infinite-dimensional intersection with the functional dimension equal to the 
difference k —p\. 

Proposition 3.12. Among multiintegrals generating the same junction stock, there 
is a unique irreducible Iq = LQ{x,y)ilj{x) (unique up to transformations Lq ^ 
Lo{x,y)R{x) with any invertible operator matrix R{x)). All other multiintegrals 
with the same function stock have the form X = Lq{x, y)S{x)(p{x) , where S{x) is an 
arbitrary operator matrix with coefficients depending on x, the number of the rows 
of this matrix is equal to its rank. 
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The proof of this Theorem repeats the proof of Proposition 12.91 and uses Theo- 
rem I3.1UI instead of Lemma 12.71 

Note that in the approach we used so far we did not allow operations with the 
arbitrary functions ifi{x). The natural admissible operations with the function sets 
ip{x) = {(pi{x), . . . ,ipkix)y are: 

• transposition of two functions (fi ^ (fj; 

• multiplication of one function with a given nonzero multiplier g{x)] 

• substitution ipj \—>- ipj + A{x)ipi, i ^ j, where A{x) is a LODO with coefficients 
depending only on x. 

These operations correspond in the representation (13. ip to: 

• transposition of the corresponding operators Li{x,y) ^ Lj{x,y) of the opera- 
tor row L{x, y) and the corresponding columns of the matrix R{x)] 

• multiplication of the corresponding operator Li{x,y) and the i-th column of 
the matrix R{x) on the right with g{x); 

• substitution Li ^ Li + LjA{x) and addition to the i-th column of R{x) of the 
j-th column multiplied on the right with A{x). 

All these operations correspond to multiplication of the both sides of the equality 
(13.11) on the right with an invertible operator matrix. It was proved in [1] that 
carrying out simultaneously the aforementioned operations with both the columns 

and the rows of the matrix R{x) lets us to reduce it to the form R{x) where Rii{x) = 
1 for 1 < i < p — 1, Rpp{x) is a LODO, and all the other elements of the matrix 
R{x) vanish. 

From this we deduce that applying to the obtained in the algorithm RMI de- 
composition L = MR the described above operations with the columns and the 
rows we get the following representation of any multiintegral: 

(Li{x,y),...,Lk{x,y)) = {Mi{x,y), . . . , Mp{x,y))R{x), 

that is Li = Ml, . . . , Lp_i = Mp_i, Lp = MpRpp and, if k > p, Lp+i = 0, . . . , 
Lk = 0. 

Thus the following Proposition is proved: 

Theorem 3.13. Using the admissible operations with the set (p{x) = {(pi{x), 
. . . ,ipk{x)y and dropping zero components, any multiintegral can be reduced to the 
form 

J = (Li{x,y), . . .,Lp_i{x,y),Mp{x,y)R{x))^{x), 
where the row (Li, . . . , Mp) gives an irreducible integral. 
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